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ABSTRACT 

Magnetic suspensions (MS) satisfy the long life 
and low loss conditions demanded by satellite and ISS 
based flywheels used for Energy Storage and Attitude 
Control (ACESE) service. This paper summarizes the 
development of a novel MS that improves reliability via 
fault tolerant operation. Specifically, flux coupling 
between poles of a homopolar magnetic bearing is 
shown to deliver desired forces even after termination 
of coil currents to a subset of “failed poles”. Linear, 
coordinate decoupled force-voltage relations are also 
maintained before and after failure by bias linearization. 
Current distribution matrices (CDM) which adjust the 
currents and fluxes following a pole set failure are 
determined for many faulted pole combinations. The 
CDM’s and the system responses are obtained utilizing 
ID magnetic circuit models with fringe and leakage 
factors derived from detailed, 3D, finite element field 
models. Reliability results are presented vs. 
detection/correction delay time and individual power 
amplifier reliability for 4, 6, and 7 pole configurations. 
Reliability is shown for two "success" criteria, i.e. (a) 
no catcher bearing contact following pole failures and 
(b) re-levitation off of the catcher bearings following 
pole failures. An advantage of the method presented 
over other redundant operation approaches is a 
significantly reduced requirement for backup hardware 
such as additional actuators or power amplifiers. 

INTRODUCTION 

Attractive magnetic bearing actuators as shown in 
Fig. 1 possess individual pole forces that vary 
quadratically with current. The net force of the bearing 
may be linearized with respect to the control voltages 
by utilizing a bias flux component. 12 Thus the X ] , *2 
and X3 forces become decoupled, i.e. dependent only 
on their respective control voltages c \^ c2^ cl) ‘ 

Maslen and Meeker 3 provided a generalization of this 
approach for heteropolar magnetic bearings (HEMB), 
which derive their bias flux from electric coils and 
utilize both N and S at different poles. 


Fault tolerant control of HEMB’s has been 
demonstrated on a 5 axis, flexible rotor test rig with 3 
CPU failures and 2 (out of 8) adjacent coil failures. 4 
Current distribution matrices for HEMB’s were 
extended to cover 5 pole failures out of 8 poles 56 and 
for the case of significant effects of material path 
reluctance and fringing. 7 

The fault tolerant approach outlined above utilizes 
a current distribution matrix (CDM) that changes the 
current in each pole after failure in order to achieve 
linearized, decoupled relations between control forces 
and control voltages, i.e. 

fcj=K vj V cj j - 1,2,3 (1) 

A failure configuration is defined by the subset of poles 
that fail due either to shorting of a turn in a coil or to 
failure of a power amplifier. In general there exist 

(2 n — 1) number of possible failure configurations for 

an n pole magnetic bearing. 

A unique contribution of the present work includes 
the extension of a CDM approach to 4, 6 and 7 pole 
homopolar magnetic bearings (HOMB). The HOMB 
commonly uses permanent magnets for its bias flux to 
increase the actuator’s efficiency and reduce heat 
generation . 8 Points on the surface of the spinning 
journal in the homopolar bearing do not experience 
north-south flux reversals thereby reducing rotor losses 
due to hysteresis and eddy currents. 

A second contribution of the present work is an 
investigation of the reliabilities of fault-tolerant 
HOMB. The reliabilities presented are system specific 
for two reasons. 

(a) An exact solution CDM may not exist for certain 
pole failure configurations. An approximate solution 
will always exist though and its effectiveness is verified 
or nullified via failure simulation for the specific 
system studied. 

(b) The two types of reliability presented correspond to 
whether a successful outcome is defined by: 

Successful outcome l (SOI): No contact between 
the shaft and catcher bearings during the failure and 
CDM implementation sequence, or 
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Successful outcome 2 (S02): Shaft contact with a 
catcher bearing then re-levitation occurs during the 
failure and CDM implementation sequence. 

Satisfaction of these success criteria will depend on the 
system studied and the delay time T d required to 

identify which poles have failed, to turn off the power 
amplifiers for these poles and to implement the 
corresponding CDM for the remaining poles. 

Two types of successful outcomes are defined in 
order to provide the system designer with magnetic 
bearing component reliabilities estimates which are 
either independent (SOI) or dependent (S02) on the 
accuracy of the catcher bearing simulation model. 
Therefore reliabilities are presented for the SOI and 

S02 conditions and for a range of f d values. 

The specific system employed for this study is a 
high speed flywheel under development for energy 
storage and attitude control applications on satellites or 
on the ISS (International Space Station). A general 
result identified from the study is an increase in 
reliability as the number of poles increase. 

FAULT TOLERANT CONTROL (FTC) 
REQUIREMENTS 

Derivation of the FTC approach requires 
applications of Ampere’s, Ohm’s, Faraday’s Laws and 
the Maxwell Stress Tensor to the multi-path magnetic 
circuit in a magnetic bearing. The physical 
requirements of FTC include 

0) Decoupling Condition: The X i control voltage 
(V ci ) does not affect the Xj control force (F x .) 

unless i = j, where Xj= X (radial) x 2 = y (radial) 
and x 3 = z (axial). 


X i 

= 0, i*j and ij = 1,2,3 (2) 

dV a 

(2) Linearity Condition: The X f control voltage ( V ci ) 
and X j control force (F xj ) are linearly related. 

F xi =K v V ci , i = 1,2,3 (3) 

where K vi is evaluated at the desired operating 
location of the shaft in the bearing. 

(3) Invariance Condition 1: The gains K v . are not 
affected by the failure. 


(4) Invariance Condition 2: The force/position gains 


dF x 

K n , ~ ~~~ L 

1 dx. 


V cJ =0,X j^X j0 ,(j=l2,3) 


/ = 1,2,3 (4) 


are not affected by the failure. The steady state 
operating point of the shaft in the bearing has 
coordinates xjq . 

The FTC requirement, Eq. (4), is automatically 
satisfied for a magnetic bearing with bias fluxes 
generated by permanent magnets (PM). This results 
since the PM’s and the resulting bias flux are 
unaffected by the failure state of the poles. 

A complete derivation of the FTC theory is 
developed next for a 6 pole homopolar combination 
(combo, radial and axial forces) magnetic bearing 
(6PHCB). The FTC theory for the 4 and 7 pole bearings 
is very similar and is not included. 

Six (6) Pole Homopolar Combo Bearing (6PHCB) 



Fig. 1 Six Pole Homopolar Combo Bearing 

Figure 1 depicts a combination (radial/axial) 6PHCB 
installed on a vertically directed shaft. The actuator has 
6 radial poles and coils and 2 axial poles and coils. The 
axial coils are wound circumferentially around the shaft 
and the radial coils are wound around the poles. The 
coil leads also form secondary coils around a common 
de-coupling choke (DC) and the axial leads also form 
tertiary coils around a second DC. The DC’s eliminate 
mutual inductances and insure that the inductance 
matrix is non-singular, which insures electric circuit 
stability. 2 The laminated construction provides for an 
accurate approximation of infinite bandwidth between 
currents and fluxes. Following common practice, the 
actuator is modeled as an equivalent circuit with 
derated magnetic strength accounting for leakage and 

derated gap flux density ( B j ) to account for fringing. 

Figure 2 shows the 6 flux paths through the radial poles 
and 2 flux paths through the axial poles. 
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Fig. 2 Equivalent Magnetic Circuit for the Six Pole 
Homopolar Combo Bearing 

The NI sources represent radial and axial control 
current flux sources. The H c L pm and R pm terms 


and N f and a ( are the number of turns on pole / and 

the gap-cross section area, respectively. The term V in 
Eq. (8) and the VI term in Eq. (7) show that the control 
flux ( VI ) varies with control current and with shaft 
position (gap values), however the bias flux ( B bias ) 

varies solely with shaft position. 

Magnetic bearings typically utilize servo power 
amplifiers (PA) that provide 1. 2-2.0 kHz bandwidth for 
inductive loads ranging between 2 mH to 8 mH. Thus it 
is acceptable to use a constant for the control current 
per control voltage gain. Let 


represent the permanent magnet source strength for 
driving bias flux,0 5 , and the reluctance of the 
permanent magnet, respectively. The magnetic circuit 
provides a useful tool to present flux conservation and 
Ampere Law relations with an equivalent electric 
circuit model. Kirchhoffs law applied to Fig. 2 yields. 



R& = NI + H (5) 

Let A represent a diagonal matrix of pole gap areas 
then by assuming uniform flux densities in each gap 


v c = (V c \ V C 2 Yc 3) T (11) 

represent the control voltages and the 8x3 matrix T is 
the current distribution matrix (CDM). Then in the 
absence of pole failures 


I' = TV C ( 12 ) 

where T includes the PA gain and the current 

distribution terms. 

Fault conditions are represented using the matrix 
K that has a null row for each faulted pole. Then the 
failed actuator control currents become 

I = KI' = KTV C (13) 

For example if coils 1 and 2 fail 

'0 0 0 0 0 0 0 o' 

00000000 /..X 

00100000 
_00010000 
”00001000 
0 0 0 0 0 10 0 

0 0 0 0 0 0 10 

0 0 0 0 0 0 0 1 


AB = < D 

(6) 

The magnetic forces are determined from the Maxwell 

B = A-'R-'N1 + A-'R-'H = VI + B bias 

(7) 

stress tensor as; 


v - a~'r~ ] n 

(8) 

Fj = B T yj B 

(15) 

B bias = A-'R-'H 

(9) 

where 


where the reluctance of gap i is 


n - diag[ ° , C0S X, / = 1 ~ 6 < n(7,7)=n(8,8) = 0 
2 Mo 

(16) 

R, = g, /(AA) 

(10) 

n =diag[ aiSin0 ‘],i = l~6’ n( 7,7) = r 2 ( 8 . 8 ) = 0 
2 Mo 

(17) 
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(7,7) = -^3 (8,8) = a ■ ■ » all other components are zero 
2 Mo 


(18) 


Substitute Eq. (13) into Eq. (7): 


B = VI + Z? 4mi = VKTVc + (19) 


where IE = FAT 

The magnetic forces are given in terms of control 
voltages and bias flux density as; 


Fj = (WVc + B bjas ) T Yj (WVc + B bka ) 

= Vc T W T y jWVc + 2B T bias VjWVc + B T bias Yj B bl 
j = 1,2,3 


( 20 ) 


The magnetic forces are proportional to the square of 
control voltages in Eq. (20). The following constraint 
equations must be satisfied in order to meet FTC 
requirements 1 (Eq. 2), 2 (Eq. 3), and 3. 


" Y \W =0 3 ^3 

2 *L^-[*vl 0 °1 

W T r 2 W = 0 3x3 

2fi L^ = [0 *v2 0] 


^^ = 03*3 


2S L^ = [o 0 at v3 ] 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 


Let 


W = [w x W 2 w 3 ] (27) 

Then the 27 constraint equations become 

wj Yj W k = 0 , /, j, k = 1,2,3 and * > / (28) 

[ 0 for i * j 

B Lr^j = \K v y/ for . = j ’V = W and j > i 

(29) 

Equations (28) and (29) are 1 8 nonlinear and 9 linear 
algebraic equations for the CDM entries, ty . The CDM 

matrix entries are obtained by requiring simultaneous 


solution of the equations in Eq. (28) and (29), and 
minimization of the Frobenius matrix norm of the CDM 
matrix. This is typically performed at the magnetic 
center, i.e. the location where the bias flux balances the 
static loads on the bearing. The norm of the current 
vector, I in Eq (13), satisfies the consistency condition 9 


IkhMMIK 

where for a Frobenius norm 


(30) 



(31) 

l‘J 


huj 2 

(32) 

iu 


£V ci 2 

(33) 


Thus by Eq. (30) reduction of |J/|| follows from 
minimizing |f|| . The Lagrange multiplier approach is 
employed to locate a solution of the equations in 
Eq. (28) and (29), that minimize ||r|| . The cost function 
is 


P 3 


1=1 7=1 k = 1 


(34) 


where p is the number of functioning poles and h k 

are the constraint equations in Eq. (28) and (29). The 
solution condition is; 

dL 


dZ m 

which implies 


— 5 ’ e {t[j * } 


(35) 


= [0 


^ 3 L 3 L 3 L 

3 /,, 3/, 7 3 / n 


3 L 3 L 3 L 

dt P i a ',3 


Of 


(36) 


The total set of equation is over-determined, i.e. more 
equations than unknowns, therefore a solution exists 
only in the least square (LS) sense. The nonlinear 
equation, LS based solver available in MATLAB is 
employed for this purpose. The effectiveness of each 
solution in satisfying the FTC requirements must be 
checked by transient response simulation of the 
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respective fault event since the LS solution is not exact. 
Fortunately the feedback control action compensates for 
the presence of residuals in the solution of Eq. (35), in 
many instances. 


These equations are solved for ty and utilizing the 

Lagrange multiplier/nonlinear least square solver 
approach discussed for the 6PHCB. 


6 Pole Homopolar Radial Bearing (6PHRB) 


Decoupling Choke 



Fig. 3 Equivalent Magnetic Circuit for the Six Pole 
Homopolar Radial Bearing 


The inductance matrix of the isolated combo 
bearing is singular because flux conservation introduces 
a dependency relation between the fluxes. This 
produces a potentially unstable operation state for the 
power amplifiers. Two decoupling chokes are added to 
the combo bearing according to the approach. 2 By 
adjusting the parameters of the decoupling chokes 

( N cl , N cl , N c2 , R cl , and R c2 ) the inductance matrix 

becomes full rank and the mutual inductances become 
zero. Similarly, a single decoupling choke is added to 
the radial bearing. 


The 6PHRB provides force solely in the two 
transverse (radial) directions. A magnetic circuit model 
for this bearing is illustrated in Fig. 3. The flux-current 
relations for this circuit are obtained by applying 
Kirchoff s laws, which yield 
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Force Linearization 

An exact solution for the ty can be obtained only 

for the no-poles failed case. Consequently the FTC 
linearization and decoupling conditions are only 
approximately satisfied and the force expressions in Eq. 
(20) are still somewhat nonlinear. Closed loop, coupled, 
flexible body simulations of the flywheel rim and shaft, 
housing, gimbals, and support structure provide 
predictions of stability, transient and steady-state 
harmonic responses. Efficient run-times for these 
models require linearized expressions for the x \ , *2 

and *3 magnetic forces. These expressions are 

obtained by applying a two-term Taylor series 
expansion about the operating point 
Po = { x j= x j0’ v cj = v cjO) • This yields 


where 


l 1 7 2 2 

R d= \Sc\d x \ x 2 

(38) 

M0 a d ’ 

The FTC requirements result in 10 

constraint 

equations 

K' r V|[(=l); t2 

(39) 

2B Ln«'=k„i o] 

(40) 

w T r2 U’=o 2x2 

(41) 

2B Lnn , ‘[» Xv 2 ) 

(42) 


- £ {“ K pij (* j */o)+ K v y {v c j v cy 0 )} (43) 
7=1 


to , , _ W?Lr,B bl J 


L =-<2 BLr,^~)k (44) 


a*,' a*, a*, 


r T m (45) 


for i=l,2,3 and j=l,2,3. The K p y and K v y 

expressions in Eq. (44) and Eq. (45) are referred to as 
’’position” and "voltage” stiffnesses respectively. The 
K v y terms are zero for i * j , only if Eq. (36) is 

satisfied exactly. Equation (20) shows that the K p y , as 

defined in Eq. (44), are independent of the ty , when 

V co is a null vector, which is typically true. 
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Flywheel and Magnetic Suspension Dynamics Model 

The novel redundant actuators operate within a 
feedback controlled system that includes both electrical 
component and structural component dynamics. A 
typical application is a flywheel module consisting of a 
high speed shaft, integrally mounted motor-generator, 
composite flywheel rim, magnetic suspension and 
flexibly mounted housing. Figure 4 depicts a module 
model with 9 rigid body structural degrees of freedom: 
rotor CG translations ( X \ r , X2 r , )* rotor rotations 

(&\ r ->&2r X housing CG translations (^ 1 / 7 ^ 2 /?) anc * 
housing rotation (# 1 / 7 , 02 / 7 )* The magnetic suspension 

employs magnetic (MB) and backup (catcher, CB) 
bearings at both the A and B ends of the module. 
Magnetic bearing clearances are approximately 0.5 mm 
so small angle motion may be assumed. Let b, d and c 
denote the magnetic actuator (Eq. 20), mass imbalance 
and catcher bearing forces, respectively. The structural 
equations of motion for the rotor are: 

M A = F A + F° + F A + F° + F A + F° i = 1.2 (46) 
M,X 3r =F£-M r g + F£ (47) 

hr #1 r +I pr °^2r = M , r + M ,,' r + M ^ (48) 
I tr 0 2r - ~ 1 pr^ir = M 2 r + M 2rc + M 2rc <49) 


For the housing the equations of motion are: 

M tf ih = Fe + Fe * ~ F ib ~ F ic ~ Fc * = 12 < 50 > 

Inhhh = M lh + M? hc + Mf hc i = l,2 (51) 

where 

M lr = ~LlF 2 A h + L B hr F 2h -L A dr F 2 A d + L B dr F 2 B d (52) 

M 2 r = 44 ~ 44 + - L B dr F B (53) 

"i * = 44 - 44 - 44 + 44 < 54 > 
M lh =-44 + 44 +i *4 -44 < 55 > 

4 = (*!/,+ Le°2h > - Q (*l A + L-e d 2h ) (56) 

4 =- K e(X\h-LeO2h)-C e {X\h-Le02h) ( 57 ) 
4 = ~X e (X2h -l40\h)-C e (X lh -L A e 6 xh ) (58) 

4 = -* e (*2A +X|^*)-C e (^2* +4%) (59) 

The simplified catcher (backup) bearing model shown 
in Fig. 5 is employed for illustrative purposes in this 
model. More sophisticated models with internal 
dynamics of races and balls or rollers are available 10 and 
could also be used in the system dynamics model. 

Let j=l,2 represent the A and B ends of the flywheel 
module in Fig. 4, respectively. Also let r represent the 
relative displacement between the catcher bearing and 
shaft at end j. 


1 / . . \2 / . . \2 


'i-fxk-x L) *[ X L- X L) 

(60) 

Fj =K c (rj -n 0 )+C c fj 

(61) 

4 = -fJ (cos Oj - n sin 6 j ) 

(62) 

4 = -Fj (sin Oj + // cos Oj ) 

(63) 
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Fig. 5 Catcher Bearing Contact Model 


Then if r 0 is the catcher bearing clearance and r ; > r$ 


Ac=H)AA 

(64) 

iL 

a 

ii 

T 

1 — 

+ 

rfc. 

(65) 


(66) 

"4- =<-■>' 

(67) 

Similarly for the axial direction if \X^ r \ > 0 


F£ =^K c {\X 2r \-r 0 )^ + C c Xi r ) 

(68) 

The mass imbalance disturbance in the model is 
described by 

A 2 

F\ d = M r ea) cos ox 

(69) 

F 2 ^ = M r ea ) 2 sin cot 

(70) 

4 = M r ea> 2 cos (ox + if/) 

(71) 


The control law utilized in the model is MIMO 
based and similar to the work of Okada" and Ahrens. 12,13 
Figure 6 illustrates the overall feedback control loop for 
the magnetic suspension. Eight power amplifiers are 
utilized for the combo bearing and 6 power amplifiers 
for the radial bearing. Five displacement sensors 
measure the relative displacements between the rotor 
and housing. Current Distribution matrices (CDM) for 
the combo and radial bearings are incorporated in the 
controllers to produce reference voltages for the 14 
power amplifiers which produce the desired currents in 
each coil. The nonlinear magnetic forces are determined 
according with Eq. (20). 


Dislurtxjnt* 



S«n*or$ 


Fig. 6 Magnetic Suspension Control Scheme 


EXAMPLES 

An example flywheel module illustrates operation 
and reliability of the redundant magnetic suspension. 
Table 1 lists the geometrical, inertia and stiffness 
parameters for the model. 

The catcher bearing contact model in Fig. 5 has a 
stiffness of 10 8 N/m, a damping of 5,000 N-s/m and a 
dynamic friction coefficient of 0.1. Table 2 shows the 
magnetic bearing parameters for the MS model. 

The ID magnetic circuit models as shown in Fig. 2 
and 3 must be adjusted to include the effects of 
recirculation leakage of the flux between the N and S 
poles of any permanent magnet and for the effect of 
non-parallel (fringing) flux flow in the air gap of each 
pole. These effects are apparent in a 3D finite element 
(FE) based simulation of the actuator as shown in 
Fig. 7. These adjustments are made with multiplicative 
factors applied to the gap flux and permanent magnetic 
(PM) coercive force in the ID model, as derived from 
the 3D FE model. The PM coercive force is derated 
from 950,000 to 514,000 in the combo bearing and 
from 950,000 to 566,000 in the radial bearing. The air 
gap flux’s are derated with a fringe factor of 0.9 for 
both the combo and radial bearings. 
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Fig. 7 3D FE Model of the Combo and Radial 6 Pole 
Actuators 

The remaining parameters of the system model include 
Displacement Sensor Sensitivity = 7874 V/m 
Displacement Sensor Bandwidth = 5000 Hz 
Power Amplifier DC Gain = 1 A/V 
Power Amplifier Bandwidth = 1200 Hz 
These 3D bearing models were also employed to verify 
the fault tolerant operation predicted with the ID 
model. An example of this is the 3 pole failure results 
shown in Table 3. The control voltage sets in this table 
are; 



\w 

0 of 

for 

set 1 

v c = (*cl Vc2 Vc-i? =' 

(0 

1 V of 

for 

set 2 


(0 

0 1 vf 

for 

set 3 


failed pole j, while the no-pole failed CDM is retained. 
The appropriate CDM for the pole-failure configuration 
being tested is then swapped in following a delay time 
x d . The MIMO control law in Fig. 6 is invariant 
throughout the entire simulation. The combo and radial 
bearing CDM’s for the no pole failed state are 


0.30789 

0.17776 

0 

0 

0.35552 

0 

-0.30789 

0.17776 

0 

-0.30789 

-0.17776 

0 

0 

-0.35552 

0 

0.30789 

-0.17776 

0 

0 

0 

-0.11530 

0 

0 

0.11530 



0.28074 

0.16209 

0 

0.32417 

-0.28074 

0.16209 

-0.28074 

-0.16209 

0 

-0.32417 

0.28074 

-0.16209 


The inductance matrix of the combo bearing with 
the two decoupling chokes is given in henries as; 


The new CDM's for the poles 1-2-3 failed case in Fig. 1 
are 



5.59x1 (T 1 

0 

0 

0 

0 

5.59x1 0" 4 

0 

0 

0 

0 

5.59x1 0’ 4 

0 

0 

0 

0 

5.59x1 0' 4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

5.59x1 (T 4 

0 

0 

0 

0 

5. 59x1 O'* 

0 

0 

0 

0 

5.85x1 0' 1 

0 

0 

0 

0 

5.83x10“ 


The inductance matrix of the radial bearing with a 
decoupling choke is given in henries as; 


3 = 


0 

0 

0 

-0.61552 

0 

0.61552 

0 

0 


0 

0 

0 

-0.37074 

-0.69571 

-0.37074 

0 

0 


0 

0 

0 

0 

0 

0 

-0.11530 

0.11530 


and 



6.76x10“* 

0 

0 

0 

0 

0 

0 

6.76x10“* 

0 

0 

0 

0 

0 

0 

6.76x10'* 

0 

0 

0 

0 

0 

0 

6.76x10“* 

0 

0 

0 

0 

0 

0 

6.76x10 * 

0 

0 

0 

0 

0 

0 

6.76x10' 


(75) 


The current produced by a power amplifier (PA) is 
turned off at the moment of failure which simulates an 
open circuit. This is implemented in the model by 
changing the K matrix in Eq. (13) from the identity 
matrix to its pole-failed value, i.e. a null row j for each 


Y D 

7 123 — 


0 

0 

0 

-0.56126 

0 

0.56124 


0 

0 

0 

-0.33768 

-0.63476 

-0.33768 


(77) 


The new CDM’s for the poles 1 -2-3-4 failed case in 
Fig. 1 are 
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rpA _ 

y 1234 ” 


o 

o 

0 

0 

-0.60742 

1.2224 

0 

0 


0 

0 

0 

0 

-1.0622 
-3.6408x1 O' 3 
0 
0 


0 

0 

0 

0 

0 

0 

-0.11530 

0.11530 


and 


rpB 

M234 — 


0 

0 

0 

0 

-0.55379 


0 

0 

0 

0 

-0.96848 


1.11460 -3.3556x10“ 


(78) 


The text below discusses two illustrative examples 
that assume identical failures in both the radial and 
combo bearings. Although this represents an extremely 
rare occurrence it serves to illustrate the method and 
analysis presented. Example 1 considers failing radial 
poles 1 , 2 and 3, and example 2 considers failing radial 
poles 1, 2, 3 and 4 in Fig. 1. Figure 8 reveals that for 
example 1 excellent control is maintained utilizing the 
no-poles failed CDM’s throughout the entire 
simulation. Consequently successful outcome criteria 
SOI is satisfied independent of the delay time x d . The 
currents in the 6 amplifiers are shown in Fig. 9 and 10 
for a failure initiation at 0.1 seconds and a delay 
duration of 20 ms. In contrast, example 2’s SOI is not 
always satisfied so that the 1 -2-3-4 poles failed CDM’s 

( 7j 234 > ^1234 ) must be activated after delay x d . The 
displacement and current responses for example 2 are 
shown in Fig. 11, 12 and 13 for x d = 20 ms. The 
reliability for example 2 will be affected by the 
selection of SOI or S02 and the delay time T d 

Successful outcome criteria 2 (S02) requires that 
the rotor successfully levitates following contact with 
the catcher bearings (CB). This is highly dependent on 
whether backward whirl (BW) develops during the 
contact period. The BW state occurs due to friction at 
the contact interface between the shaft and CB, which 
forces the shaft to whirl (precess) in a direction 
opposite to the spin direction. Figure 14 shows an 
example of this state with p=0.3, C c =10 5 N-s/m and 
K=10 8 N/m. The BW eccentricity is the CB clearance 
(typically 0.25 mm) for a rigid rotor, and possibly a 


much larger value for a flexible shaft. The whirl 
frequency typically ranges from 0.4- 1.0 times the spin 
frequency. This creates a potentially large centrifugal 
force that can damage the CB’s or deflect the shaft into 
the MB’s. The BW condition is mitigated by proper 
design of the flexible damped support, preload, 
clearance and friction coefficient for the CB’s. Re- 
levitation off of the CB’s is very difficult once BW has 
fully developed. Figure 15 shows the displacements for 
a successful re-levitation event with poles 1 -2-3-4 
failed, x d =100 ms and C =5000 N-s/m. 



0 0.05 0.1 0.15 0.2 

Time (s) 


Fig. 8 Rotor Displacements in the Radial and Axial 
Direction for Example 1 



Tim* <*) Tim* W) 

Fig. 9 Current Responses in Combo Bearing for 
Example 1 
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Axial X3 (mm) Radial (mm) Racial X1 A (mm) 


Currents in Radial Bearing 



0 0.05 0.1 0.15 0.2 0 0.05 0.1 0.15 0.2 



0 0.05 0.1 0.15 0.2 0 0.05 0.1 0.15 0.2 



Time (s) Time (s) 


Fig. 10 Current Responses in Radial Bearing for 
Example 1 



0 0.05 0.1 0.15 0.2 


Time (s) 

Fig. 1 1 Rotor Displacements in the Radial and Axial 
Direction for Example 2 



*n»|l] 


Fig. 12 Current Responses in Combo Bearing for 
Example 2 


Currents in Radial Bearing 



0 0.05 0.1 0.15 0.2 0 0.05 0.1 0.15 0.2 



0 0.05 0.1 0.15 0.2 0 0.05 0.1 0.15 0.2 


Time (s) Time (s) 

Fig. 13 Current Responses in Radial Bearing for 
Example 2 



-0,3 -0.2 -0.1 0 0.1 0,2 0.3 


XI A (mm) 

Fig. 14 Orbit Plot of the rotor at C.B. 1 with 
C c — 10 5 Nt-m/s 



0 0.1 0.2 0.3 

Tim® (s) 


| 0.2 
ro n 

* \ 

* -0 2 

0 0.1 0.2 0.3 

Time (s) 

Fig. 15 Rotor Displacements in the Radial and Axial 
Direction During Successful Re-levitation 
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The reliability of a magnetic suspension (MS) is 
determined by considering the number of failed pole 
states that still meet the SOI or S02 criteria. This is 
dependent on the time delay x d , modeling assumptions, 
number of poles in the bearing and the reliability of the 
power amplifier/coil units that drive and conduct the 
bearing currents. The 4 pole and 7 pole configurations 
require 2 less or 1 more power amplifiers than the 6 
pole configuration, respectively. The radial pole and 
permanent magnet cross-section areas, the number of 
turns of each radial coil, and the coercive force and the 
length of the permanent magnets for the 4 and 7 pole 
bearings are identical to those of the 6 pole bearing. 

The no-pole failed CDM’s for the 7 pole bearing 
are 


’ 0.33071 

-0.017910 

-4.6780x1 0' 3 

0.17799 

0.26601 

3.6330x1 0~ 3 

-0.067428 

0.26396 

-2.7850x1 0' 3 

-0.26676 

0.16181 

7.8561xl(T 4 

-0.29880 

-0.15575 

1.7079x1 0" 3 

-0.038458 

-0.29204 

-4.2590x1 0' 3 

0.15965 

-0.23855 

4.959 lxl 0" 3 

0 

0 

-0.099552 

0 

0 

0.099552 


0.28402 

0 

0.17708 

0.22206 

-0.063201 

0.27690 

-0.25589 

0.12323 

-0.25589 

-0.12323 

-0.063201 

-0.27690 

0.17708 

-0.22206 


The no-pole failed CDM’s for the 4 pole bearing are 


' 0.52550 

0 

0 

0 

0.52550 

0 

-0.52550 

0 

0 

0 

-0.52550 

0 

0 

0 

-0.17043 

0 

0 

0.17043 


0.46539 
7.3028x1 O' 4 
-0.46403 
7.3028x1 0~ 4 


-1.07319xl0" 3 
0.46371 
-1.073 19x1 0” 3 
-0.46571 


and 


(80) 


The radial pole failure simulations are conducted 
with the combo bearing operating in a no-pole failed 
state, and vice versa. Failure occurs at 0.1 seconds into 
the simulation and swapping in of the new CDM occurs 
at a delay time x d later. The number of j unfailed pole 
cases for an n pole bearing is given by the formula 


/ 


»j 


L_ 

J) Mn~j)\ 


(81) 


Table 4 summarizes the results of these simulations for 
swapping in the appropriate poles-failed (new) CDM 
for the delay times x d of 20, 60 and 100 ms, 
respectively. The S01+S02 column considers all cases 
when either SOI or S02 occurs. 

An n pole, fail-safe, homopolar magnetic bearing is 
similar to a m-out-of-n system in a reliability model if 
stable control is maintained (SOI or S02) when at 
minimum m of the n poles (P.A. plus coil) are unfailed. 
Let R p represent the reliability of a "pole”, i.e. of the 
power amplifier plus its pole coil, at some specific point 
in its expected lifetime. Also assume that "poles" are 
identical and act independently. The system reliability 
then become 


R s = i a k R p k {\-R p ) n ~ k (82) 

k-m 

where (X k are the number of SOI (or S01+S02) cases 

in Table 4. The integer m in Eq. (82) is the minimum 
number of unfailed poles that are required for the n pole 
bearing to successfully levitate the shaft. The (n,m) 
pairs determined in this example are (4,2), (6,2) and 
(7,2). Figures 16-18 show system reliability vs. R p plots 
for the 4, 6 and 7 pole radial bearings for SOI and 
(S01+S02) and x d equal to 20, 60 and 100 ms, 
respectively. Figures 19-21 repeat Fig. 16-18 for the 
zoomed-in range 0.9<R p <l. Similarly, Fig. 22-24 show 
the system reliabilities for radial actuation of the combo 
bearings. Figure 25-27 repeat Fig. 22-24 for the 
zoomed-in range 0.9<R p <l. Axial control reliability is 
not considered in these figures since it is typically 
independent of radial direction control. 

Table 5 and 6 show system reliabilities R s for 
several values of pole reliability R p . The delay time x d , 
number of poles n, and success criterion are also varied 
in these tables. A particular notable result shown in 
these tables is the high reliability of the 4, 6 and 7 pole 
bearings even if the CDM is not changed after failure. 
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Fig. 16 System Reliabilities of 4, 6 and 7 Pole Radial Fig. 19 System Reliabilities of 4, 6 and 7 Pole Radial 
Bearings with x d = 20 ms Bearings with x d = 20 ms 



Fig. 17 System Reliabilities of 4, 6 and 7 Pole Radial 
Bearings with x d = 60 ms 



Fig. 18 System Reliabilities of 4, 6 and 7 Pole Radial 
Bearings with x d = 100 ms 



Fig. 20 System Reliabilities of 4, 6 and 7 Pole Radial 
Bearings with x d = 60 ms 



Fig. 21 System Reliabilities of 4, 6 and 7 Pole Radial 
Bearings with x d = 100 ms 
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Fig. 22 System Reliabilities of 4, 6 and 7 Pole Combo Fig. 25 System Reliabilities of 4, 6 and 7 Pole Combo 
Bearings with x d = 20 ms Bearings with x d = 20 ms 



Fig. 23 System Reliabilities of 4, 6 and 7 Pole Combo Fig. 26 System Reliabilities of 4, 6 and 7 Pole Combo 
Bearings with x d = 60 ms Bearings with x d = 60 ms 



Fig. 24 System Reliabilities of 4, 6 and 7 Pole Combo Fig. 27 System Reliabilities of 4, 6 and 7 Pole Combo 
Bearings with x d = 100 ms Bearings with x d = 100 ms 
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Table 1 Flywheel Model Parameters 


M r 

29.644 kg 

M h 

34.428 kg 

hr 

0.26233 kg-m 2 

V 

0.11129 kg-m 2 

h \h 

1.5337 kg-m 2 

hih 

1.3993 kg-m 2 

K e 

3.5024x10 s N/m 

c* 

5.2535 x10 s kg/s 

0) 

60,000 rpm 

e 

2.54X10 -6 m 

r A 

L br 

0.14051 m 

l b 

L br 

0.13360 m 

jj & 
L dr 

0.14051 m 

l b 

L dr 

0.13360 m 

[A 

^ sr 

0.17846 m 

l b 

.ST 

0.16974 m 

[A 

^ cr 

0.26765 m 

I? 

^ cr 

0.28067 m 

l a 

L bh 

0.14051 m 

i B 

L bh 

0.13360 m 

l a 

L sh 

0.17856 m 

I? 

L sh 

0.16974 m 

l a 

ch 

0.26765 m 

l b 

0.28067 m 

l a 

0.26765 m 

l b 

0.28067 m 

V 

n 12 




Table 2 Magnetic Bearing Parameter List 



Combo Bearing 

Radial Bearing 

Air Gap 

radial: 5.080 XlO^/W 

radial: 5.080 X KT^m) 

axial: 5.080xl0^tf7 

dead pole: 0.00203(/w) 

Radial Pole Face Area 

3.924x10' 4 (w 2 ) 

4.746x10" 4 (/» 2 ) 

Axial Pole Face Area 

1.719xl0" 3 (m 2 ) 

N/A 

Dead Pole Face Area 

N/A 

4.962xl0" 3 (m 2 ) 

Total Face Area of PM 

3.178x10' 3 (w 2 ) 

3.844x10" 3 (/w 2 ) 

Length of PM 

0.0101(m) 

O.OlOl(m) 

Number of Turns of Radial Coil 

24 

24 

Number of Turns of Axial Coil 

37 

N/A 

Relative Permeability of PM 

1.055 

1.055 

Coercive Force of PM 

950000 (A/m) 

950000 (A/m) 
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Table 3 ID and 3D Model Comparison of Predicted Forces for 6 Pole Combo Bearing 


Control 

Voltage 

Set 

Force 

Direction 

Force (N) 

No Poles Failed 

3 Poles Failed (Poles 1-2-3 

ID Model 

3D Model 

ID Model 

3D Model 

1 

XI 

11.64 

12.95 

11.64 

12.96 

1 

X2 

0 

0.01 

-0.14 

-0.25 

1 

X3 





2 

XI 





2 

X2 

11.64 

13.3 

11.59 

13.17 

2 

X3 

0 

0.08 

0 

-0.05 

3 

XI 

0 

-0.4 

0 

-0.4 

3 

X2 

0 

0.66 

0 

0.66 

3 

X3 

8.9 

9.4 

8.9 

9.4 


Table 4 Summary of Simulation for Reliability Study 


n pole 
Bearing 

Failed 

Bearing 

No. of 
unfailed 
Poles (j) 

No. of 
Simulation 

(Inj) 

Delay time T d 
20 ms 

Delay time T d 
60 ms 

Delay time 
100 ms 

No. of 
SOI 
cases 

No. of 
SOl+SO 
2 cases 

No. of 
SOI 
cases 

No. of 
SOl+SO 
2 cases 

No. of 
SOI 
cases 

No. of 
SOl+SO 
2 cases 

4 

Radial 

2 

6 

4 

4 

4 

4 

4 

4 

3 

4 

4 

4 

4 

4 

4 

4 

4 

1 

i 

i 

1 

1 

i 

i 

Combo 

2 

6 

4 

4" 

4 

4 

4 

4 

3 

4 

4 

4 

4 

4 

4 

4 

4 

1 

1 

1 

1 

1 

1 

1 

6 

Radial 

2 

15 

12 

12 

6 

12 

0 

12 

3 

20 

20 

20 

9 

20 

8 

20 

4 

15 

15 

15 

15 

15 

15 

15 

5 

6 

6 

6 

6 

6 

6 

6 

6 

1 

1 

1 

1 

1 

1 

1 

Combo 

2 

15 

12 

12 

1 

12 

0 

11 

3 

20 

20 

20 

15 

20 

9 

20 

4 

15 

15 

15 

15 

15 

15 

15 

5 

6 

6 

6 

6 

6 

6 

6 

6 

1 

1 

1 

1 

1 

1 

1 

7 

Radial 

2 

21 

16 

21 

0 

21 

0 

21 

3 

35 

33 

34 

11 

29 

4 

30 

4 

35 

35 

35 

29 

35 

20 

35 

5 

21 

21 

21 

21 

21 

21 

21 

6 

7 

7 

7 

7 

7 

7 

7 

7 

i 

1 

1 

1 

1 

1 

1 

Combo 

2 

21 

13 

14 

0 

14 

0 

14 

3 

35 

28 

28 

13 

28 

5 

28 

4 

35 

35 

35 

28 

35 

21 

35 

5 

21 

21 

21 

21 

21 

21 

21 

6 

7 

7 

7 

7 

7 

7 

7 

7 

i 

1 

1 

1 

1 

1 

i 
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Table 5 System Reliabilities of Radial Bearings vs. R p , n, t, and Success Criterion 


R 

p 

ms 

SOI or 
S01+S02 

R for CDM Swap 

ped 

F R for No-Poles Failed CDM 

n=4 

n=6 

n=7 

n=4 

n=6 

n=7 

0.900 

20 

SOI 

0.98010000 

0.99970200 

0.99980730 

0.98010000 

0.98998200 

0.98772210 

60 

SOI 

0.98010000 

0.99071100 

0.99413730 

\ 100 

SOI 

0.98010000 

0.98998200 

0.98772210 

20 

S01+S02 

0.98010000 

6.99970200 

0.99992070 

60 

S01+S02 

0.98010000 

0.99970200 

0.99955620 

100 

S01+S02 

0.98010000 

0.99970200 

0.99962910 

0.950 

20 

SOI 

0.99500625 

0.99998128 

0.99998777 

0.99500625 

0.99862753 

0.99830066 

60 

SOI 

0.99500625 

0.99873470 

0.99925449 

100 

SOI 

0.99500625 

0.99862753 

0.99830066 

20 

S01+S02 

0.99500625 

0.99998128 

0.99999454 

60 

S01+S02 

0.99500625 

0.99998128 

0.99996774 

100 

S01+S02 

0.99500625 

0.99998128 

0.99997310 

0.990 

20 

SOI 

0.99980001 

0.99999997 

0.99999998 

0.99980001 

0.99998821 

0.99998529 

60 

SOI 

0.99980001 

0.99998918 

0.99999400 

100 

SOI 

0.99980001 

0.99998821 

0.99998529 

20 

S01+S02 

0.99980001 

0.99999997 

0.99999999 

60 

S01+S02 

0.99980001 

0.99999997 

0.99999994 

100 

S01+S02 

0.99980001 

0.99999997 

0.99999995 


Table 6 System Reliabilities for Radial Actuation of the Combo Bearings vs. R p , n, t d and Success Criterion 


R , 

ms 

SOI or 
S01+S02 

R„ for CDM Swapped 

I R. for No-Poles Failed CDM 

n=4 

n=6 

n=7 

n=4 

n=6 

n=7 

0.900 

20 

SOI 

■imiitiiniil 

0.99970200 

0.99941850 

0.98010000 

0.990711 

0.98845110 

60 

SOI 


0.99516600 

0.99362700 

100 


0.98010000 

0.99071100 

0.98845110 

20 


0.98010000 

0.99970200 

0.99942660 

60 

ESEEEBE1 

0.98010000 

0.99970200 

0.99942660 

r 100 

mSBESM 

0.98010000 

0.99962100 

0.99942660 

0.950 

20 

SOI 


■rami** 

0.99996013 

0.99500625 

0.99873470 

0. 99840783 

60 

SOI 

— asa 


0.99916339 

100 

SOI 

0.99500625 

0.99873470 

0.99840783 

20 

WBBESM 

0.99500625 

0.99998128 

0.99996041 

60 


0.99500625 

0.99998128 

0.99996041 

100 

S01+S02 

0.99500625 

0.99997564 

0.99996041 

0.990 

20 

SOI 

0.99980001 

0.99999997 

0.99999993 

0.99980001 

0.99998918 

0.99998626 

60 

SOI 

0.99980001 

0.99999501 

0.99999306 

100 

SOI 

0.99980001 

0.99998918 

0.99998626 

20 

S01+S02 

0.99980001 

0.99999997 

0.99999993 

60 

S01+S02 

0.99980001 

0.99999997 

0.99999993 

100 

S01+S02 

0.99980001 

0.99999996 

0.99999993 
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SUMMARY AND CONCLUSION 

This manuscript presents a description, analysis, 
example and reliability study for novel, redundant, 
radial and combination, homopolar magnetic bearings. 
Current distribution matrices (CDM) are evaluated 
based on the set of poles that have failed and the 
requirements for uncoupled force/voltage control, 
linearity and specified force/voltage gains that are 
unaffected by the failure. The CDM algorithm also 
determines the CDM with a minimum Frobenius norm 
which provides reduced effort (current required) 
operation of the HOMB. An advantage of the HOMB 
over a HEMB is the automatic invariance of the 
position stiffness before and after pole failure. This 
results from the bias flux source being permanent 
magnets. A simplified catcher bearing model is 
employed to evaluate the improvement in reliability 
which results from utilizing a success criterion (S02) 
based on re-levitation after catcher bearing contact vs. a 
criterion (SOI) which excludes all contacts with the 
spinning shaft. The SOI criterion is more conservative 
since it does not depend on the accuracy of the catcher 
bearing model used in the simulation. 

The numerical example presented exhibits several 
interesting trends which include 

the reliability of the 4, 6 or 7 pole bearing is high 
even if the reliability of the pole decreases with 
time to 0.90. 

increased reliability with increased number of 
poles 

high reliability without replacing the no-poles 
failed CDM with the appropriate poles-failed CDM 
successful levitation with only 2 unfailed poles for 
the n = 4, 6 and 7 pole HOMB's 
successful fault tolerant operation without changes 
to the MIMO control in Fig. 18. 

Future work includes tests of 6 pole radial and 
combination HOMB's in the G2-ACESE flywheel 
module at NASA Glenn. A higher fidelity catcher 
bearing model will also be employed in future 
simulations to provide a more realistic evaluation of the 
S02 type reliabilities. Finally the MIMO control will 
also be modified to improve the reliability without 
replacing the no-poles failed CDM with a pole-failed 
CDM's. 

NOMENCLATURE 

a = radial pole face area of combo 

bearing 

f 

a = axial pole face area of combo 

bearing 


a d 

— 

dead pole face area of radial 
bearing 

B 

= 

flux density vector 

F bias 

= 

bias flux density vector 

C c 

= 

contact damping 

C e 

= 

housing damping 

e 

= 

rotor eccentricity 

F 0' F 0 

= 

force 

g 

= 

acceleration of gravity 

S() 

= 

air gap of the radial pole 

Sod 

H c 

— 

air gap of dead pole of radial 
bearing 

‘ 

coercive force of permanent 
magnet 

h) 


current 

J >( )’ 7 p( ) 

= 

transverse and polar moment of 
inertia 

*,(X) 

= 

position stiffness 

K v { )( ) 

= 

control voltage stiffness 

K c 

= 

contact stiffness 

K e 

= 

housing stiffness 

L nm 

pm 

= 

length of permanent magnet 

tA tB 

H )'H ) 

= 

distance from the center of 

M r 


flywheel or housing coordinate 

= 

mass of rotor 

M h 

= 

mass of housing 

<)•<) 

= 

moment 

*<> 

= 

number of turns of coils 

N cl ,N c2 ,N ci 

*, 

= 

number of turns on decoupling 
chokes 

= 

reluctance of air gap 

R pm 

= 

reluctance of permanent magnet 

RcX ’ F c2 

= 

air gap reluctance of decoupling 
chokes 

To Jo 

= 

current distribution matrix 
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= flux vector 

= rotor angular velocity 
= dynamic friction coefficient 
= angle of the radial pole 
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